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FOREWORD 


This  report  is  the  first  of  a  series  dealing  with  a  theoretical 
study  of  the  phenomenology  of  bumper  protection  from  meteoric  attack.  It 
was  prepared  by  Dr.  G-  V.  Bull  of  McGill  University,  Montreal,.  Quebec, 
serving  as  a  consultant  to  Arthur  D.  Little,  Inc.,  in  connection  with 
NASA  Contract  No.  NAS  5-664. 

Dr.  Bull  is  responsible  for  the  approach  to  the  problem  and  the^ 
draft  of  this  report.  It  is  inevitable  that  various  aspects  will  emerge  in 
altered  emphasis  subsequent  to  technical  discussion  with  others.  In  arrang¬ 
ing  for  the  Issuance  of  this  report  from  Arthur  D.  Little,  Inc.,  it ^seemed 
appropriate  that  one  should  attempt  to  update  this  work  so  as  to  reflect 
those  considerations  which  emerged  from  technical  discussions  among  Dr. 

Bull,  Dr.  Sydney  Goldstein  of  Harvard  University,  and  the  writer.^  The  re¬ 
port  as  prepared  by  Dr.  Bull  is  reproduced  with  but  minor  alterations;  the 
various  fe-iures  of  our  discussions  which  it  appeared  appropriate  to  record 
are  presented  as  notes  in  this  foreword.  In  the  interest  of  brevity  the 
various  definitions  of  symbols,  etc.,  have  not  been  restated  in  the  fore¬ 
word.  For  this  reason  it  is  suggested  that  one  read  the  main  text  and  re¬ 
fer  to  these  notes  where  such  references  are  suggested  in  the  text.  Sup¬ 
plemental  passages  added  to  the  body  of  the  main  text  in  the  course  of 
preparation  for  its  issuance  as  an  Arthur  D.  Little,  Inc.,  report  appear 
in  double  parenth  ses. 

Note  !•  Most  of  the  materials  of  interest  will  probably  exhibit  a  poly- 
tropic' constant  in  excess  of  2.  In  this  case  the  expansion  occasioned  by 
the  rarefaction  wave  will  be  somewhat  constrained  by  the  bumper  since  the 
head  of  the  rarefaction  wave  will  not  emerge  behind  the  bumper  plate  (this 
is  discussed  in  the  text).  Attention  must  also  be  given  to  the  possi¬ 
bility  that  radial  decompression  of  the  pellet  may  occur  in  front  of  the 
bumper  plate.  This  circumstance  would  require  that 


7p  +  1 


1) 


(as  is  readily  shown  on  combining  equations  l4,  15,  and  l6  of  the  text). 

Such  a  circumstance  might  well  be  important  if  Pp^/pbo  «  current 

meteoric  e\ldence  would  suggest. 

Note  2:  Chronologically  Section  VII  was  completed  prior  to  a  substantial 

reappraisal  of  Section  VI.  As  a  consequence,  the  model  on  which  the  expan¬ 
sion  core  is  calculated  is  appropriate  only  for  7p  <  2.  Values  for  7b  >  ^ 
are  shown  as  dashes,  and  one  would  expect  the  core  angle  to  be  overestimated 
by  this  process.  The  assumption  that  the  appropriate  7  for  the  expansion 
process  is  the  same  as  that  for  the  compressive  phase  might  be  questioned; 
but  such  refinements  must  await  subsequent  study.  R.  H.  J. 


Arthur  ZB.lCittlc.Knr. 
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ABSTRACT 


This  note  outlines  the  results  of  a  short  theoretical  investiga¬ 
tion  into  the  impact  of  pellets  with  thin  plates  at  very  high  velocities. 

In  order  to  render  the  problem  tractable  to  analysis,  a  model  has  been  as¬ 
sumed  which  allows  use  of  essentially  one -dimensional  theory  to  predict 
both  the  impact  shock-induced  conditions  and  the  subsequent  expansion  flow, 
although  the  latter  state  has  been  treated  subsequently  as  a  general  axi- 
symmetric  three  dimensional  flow.  The  problem  has  been  considered  in 
three  phases. 

The  first  phase  is  the  shock-induced  impact  states,  which  for 
thin  bumper  plates  permits  one-dimensional  theory  to  be  reasonably  appli¬ 
cable  .  Based  on  the  concept  of  the  strong  shock -wave,  it  is  seen  that 
the  initial  impact -induced  states  may  be  expressed  readily  in  terms  of  im¬ 
pact  velocities,  initial  density  ratios,  and  the  polytropic  gas  exponent 
in  the  shocked  states. 

The  second  phase  as  discussed  in  this  note  results  from  the  ex¬ 
tremely  high  temperatures  generated  on  impact .  The  resultant  radiative 
loss  mechanisms  result  in  considerable  reduction  in  the  thermodynamic 
states  from  which  expansion  occurs.  Detailed  treatment  of  this  phase 
will  appear  in  subsequent  notes. 

The  third  phase,  the  break  up  of  the  pellet  into  an  expanding 
particle  cloud,  has  been  approached  in  this  note  from  the  restricted 
point  of  view  of  one-dimensional  shock  and  rarefaction  waves  interacting 
with  a  cylindrical  rarefaction.  This  non-stationary  wave  system  tends  to 
accelerate  the  outer-material,  essentially  redistributing  the  initial 
kinetic  energy  of  the  pellet.  Only  brief  consideration  into  the  bounda¬ 
ries  of  this  flow  has  been  made  in  the  present  note.  Subsequent  detailed 
treatment  of  the  expanding  flow  phase  will  be  presented  in  the  near  fu¬ 
ture  .* 


See  Foreword, 


Note  1 


Arthur  21.H.ittle,Knr. 
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I .  INTRODUCTION 


Experimentally,  it  has  been  observed  that  the  impact  of  small 
bodies  with  thin  plates  at  hypervelocity  speeds  results  in  the  complete 
disintegration  of  the  pellet  after  passage  through  the  plate .  An  expand¬ 
ing  cloud  of  particles  is  obtained,  the  shape  of  this  cloud  being  ex¬ 
tremely  sensitive  to  the  geometry  of  the  impact  as  illustrated  in  Figures 
1  and  2  where  an  end-on  impact  of  a  cylindrical  pellet  results  in  a  more 
or  less  circular  cross  section  of  pellet  spray  while  the  same  model  col¬ 
liding  at  approximately  ^5°  results  in  the  mushroom-type  cross  section  shown 
in  Figure  2.  The  nature  of  the  particle  cloud  is  among  other  things,  a  func¬ 
tion  of  impact  velocity;  at  low  enough  velocities  the  pellet  will  be  frag¬ 
mented  with  some  pulverization,  whereas  at  somewhat  higher  velocities,  com¬ 
plete  pulverization  may  occur,  and  finally  at  extremely  high  velocities, 
complete  vaporization  of  the  pellet  and  displaced  bumper  material  may  be 
expected.  The  little  experimental  evidence  available  (all  at  relatively 
low  velocity)  tends  to  indicate  that  the  spray  cross  section  is  not  mark¬ 
edly  dependent  upon  the  impact  velocity. 

The  impact  of  small  pellets  into  essentially  semi-infinite  tar¬ 
gets  has  been  the  subject  of  study  for  many  years;  the  extreme  complexity 
of  the  theoretical  problem  has  tended  towards  reliance  on  a  general  semi- 
empirical  approach  (References  1,2).  The  theoretical  approaches  have  been 
based  upon  a  fluid  model,  since  interest  in  general  hinges  around  the  case 
when  shock  pressures  are  at  least  an  order  of  magnitude  greater  than  yield 
stresses.  Bjork  (Reference  3)  has  numerically  solved  the  hydrodynamic 
equations  to  determine  crater  volu-me  dependence  on  impact  velocity. 
Stanyukovich  (Reference  4),  in  considering  the  impact  of  large  meteors  on 
planets  having  no  atmospheres,  has  noted  the  strong  similarity  to  explo¬ 
sions  providing  the  collision  kinetic  energy  is  greater  than  the  energy 
density  of  the  crystal  lattice.  Even  when  the  collision  energy  is  less 
than  required  for  vaporization,  he  has  noted  that  pulverization  can  occur, 
resulting  in  a  quasi-fluid  flow.  The  \^porization  energy  for  most  materi¬ 
als  (sand,  clay,  granite,  aluminum,  steel,  etc.)  is  of  the  order  of  10^® 
to  10^2  ergs  per  gram.  For  velocities  in  excess  of  10  km  per  second,  va¬ 
porization  would  occur  and  the  conditions  behind  the  impact  shocks  are 
similar  to  those  behind  explosive  shocks. 

In  any  model,  the  displacement  phenomena  are  initiated  by  shocks 
generated  at  the  instant  of  impact;  these  shocks  are  propagated  through 
the  collided  and  colliding  bodies.  In  a  series  of  experimental  studies  on 
the  equation  of  state  conducted  at  the  Los  Alamos  Scientific  Laboratory, 

Walsh  et  al  (References  5>6,7)  used  a  large  plane-wave  explosive  accelera¬ 
tor  to  launch  a  0;22  cm  thick  steel  or  nickel  driver  plate  against  desired 
specimens  in  the  5  to  5  km/sec  velocity  range.  The  subsequent  shock  wave 
and  free -surface  velocity  was  measured,  and  from  the  mechanical  conserva¬ 
tion  conditions,  pressure  and  density  in  the  shocked  state  were  determined. 

In  these  studies,  shock  pressures  of  up  to  two  megabars  were  determined. 

While  this  is  below  the  pressures  experienced  in  meteor  impact,  it  does 
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FIGURE  1  Sketch  of  Spray  Pattern  from  End-on  Collision  with 

.  020"  Al  Bumper  Plate  (!?'  ahead  of  witness  plate)  •  f 

Zelux  Plastic  Cylinder  .550"  Dia. ,  350"  Long, 

6  km/sec.  Semi  Vertex  Spray  Angle  20°  | 
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FIGURE  2  Photograph 

Spray  Pattern  from  a  0.55"  Dia.  Plastic  Pellet  (Cylinder), 
Striking  at  Approximately  45“ 

.020"  Al  Bumper  Plate  12"  ahead  of  Witness  Plate,  Cylinder 
Length  0.350" 

6  Kilometers/sec. 
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provide  reliable  experimental  data  on  shock  propagation  in  solids  in  circimi- 
stances  such  that  vaporization  is  negligible.  The  problem  of  practical  in¬ 
terest  here  concerns  the  determination  of  protection  effectiveness  of  thin 
bumper  plates  against  meteor  impact.  Pictorially  speaking,  we  may  imagine 
that  subsequent  to  iii5)act,  the  colliding  pellet  and  entrained  bumper  mate¬ 
rial  explode,  imparting  lateral  momentum  and  energy  to  particles  of  the 
pellet  at  the  expense  of  the  original  uni -directional  pellet  momentum.  The 
spreading  of  the  momentimi  and  energy  over  a  relatively  large  area,  combined 
with  the  transfer  into  lateral  motion,  renders  the  pellet  relatively  incap¬ 
able  of  further  penetration,  the  exact  nature  of  the  damage  inflicted  on 
any  body  behind  the  bumper  being  dependent,  among  other  things,  on  the  sepa¬ 
ration.  The  central  theoretical  problem  to  be  considered  in  this  and  sub¬ 
sequent  reports,  is  the  determination  of  the  translational  energy  distribu¬ 
tion  in  the  spray  cross  section  as  a  function  of  the  pertinent  parameters, 
viz:  distance  back  of  bumper,  impact  velocity,  pellet  and  bumper  material, 
etc . 


The  attack  on  this  problem  will  consist  first  in  determining  the 
initial  impact  induced  states  on  the  basis  of  one -dimens ional  theory,  (end- 
on  impact  of  a  simple  cylinder).  Explicit,  tractable,  analytical  expres¬ 
sions  will  be  sought  to  express  the  shocked  states  in  tenns  of  the  appro¬ 
priate  initial  parameters  (biimper/pellet  density  ratio,  specific  heat  ratio, 
and  impact  velocity).  Radiative  loss  processes  will  be  delineated  to  es¬ 
tablish  the  dominant  mechanism.  Subsequently,  the  development  of  a  three- 
dimensional  expansion  (explosion-like)  profile  based  on  the  propagation 
of  appropriate  rarefaction  waves  originating  at  the  free -boundaries  will 
be  considered  on  the  basis  of  self-similar  motions.  Finally,  the  treat¬ 
ment  will  be  extended  to  other  shapes,  and  porous  materials.  This  work 
has  been  iindertaken  as  part  of  the  Arthur  D.  Little,  Inc.,  program  on 
space  storage  of  propellants,  a  program  being  conducted  for  NASA.  From 
time  to  time  reports  will  be  issued  covering  completed  phases  of  the  work; 
this  is  the  first  such  note  and  has  as  its  purpose  the  defining  of  the 
problem  and  method  of  attack. 
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II.  THE  IMPACT  MODEL 


The  impact  model  may  properly  he  considered  in  two  distinct 
phases.  The  first  phase  establishes  the  initial  shock  mechanisms  and 
states  along  the  lines  of  conventional  (References  ^,8)  one  dimensional 
hydrodynamic  shock  propagation.  To  render  the  problem  tractable,  the  as¬ 
sumption  of  strong  shock  wave  theory  is  made  in  this  work,  that  is,  the 
density  ratio  is  assumed  to  have  reached  a  limiting  value.  While  this 
cannot  be  strictly  true,  it  should  provide  a  reasonable  approximation, 
and  enables  the  derivation  of  simple  expressions  for  the  thermodynamic 
states  generated  on  impact.  The  second  phase  follows  the  propagation  of 
these  initial  planar  fronts  and  the  explosion-like  expansion  flow  of  the 
vaporized  pellet  material. 

The  impact  model  assumed  is  illustrated  in  Figure  3-  Condi¬ 
tions  just  at  impact  are  illustrated  in  sketch  (a);  the  cylinder  of  ra¬ 
dius  Tq  strikes  the  bumper  plate  at  velocity  V,  the  bumper  plate  thick¬ 
ness  6  being  much  less  than  rp.  The  impact  results  in  the  non-stat ionary 
wave  system  illustrated  in  (b);  a  transmitted  shock  Sg  proceeds  into  the 
bumper  plate  at  velocity  wg,  while  a  reflected  shock  proceeds  into  the 
pellet  at  v’elocity  wp.  The  physical  functions  of  these  shocks  are  such 
as  to  establish  the  contact  interface  in  the  ideal  mechanical  model  as¬ 
sumed  (i.e.  that  there  is  no  penetration  mixing  of  bumper  and  pellet  mate¬ 
rial).  Thus  at  the  contact  plane  the  bumper  shock  accelerates  the  bumper 
material  to  a  velocity  equal  to  the  decreased  pellet  velocity  resulting 
from  gp.  This  is  defined  analytically  by  the  relation: 


V  -  Up  =  u^ 


(1) 


((This  simply  implies  that  the  pellet  and  bumper  remain  in  contact.)) 
Also  we  have 


P 


B 


(2) 


((I.e.  there  is  no  discontinuity  of  pressure  across  the  contact  zone.)) 
The  contact  surface  is  then  separating  bumper  and  pellet  material  and 
traveling  to  the  right  at  the  particle  velocity.  In  many  instances  (as 
vill  emerge  from  the  analysis)  all  these  waves  will  travel  at  velocities 
considerably  less  than  the  impact  velocity  V,  in  which  case  all  will  be 
swept  to  the  right  with  respect  to  an  axis  system  fixed  in  the  bumper. 
((In  Section  VI,  Pellet  Dynamics,  it  is  pointed  out  that  the  rarefaction 
wave  in  the  bumper  will  move  to  the  left  with  respect  to  an  axis  system 
fixed  in  the  bumper,  provided  73  exceeds  2.  Note  1  of  the  foreword  indi¬ 
cates  the  circumstance  under  which  the  compression  of  the  pellet  may  move 
to  the  left  relative  to  axes  fixed  in  the  bumper.)) 
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(a) 

CONDITIONS  AT  IMPACT 


REFLECTION  OF  BUMPER  SHOCK 
AS  A  RAREFACTION.  B 


TRAVELS  AT  VELRCITV  Rp  WITH  RESPECT 
TO  PELLET,  *p  <  V.  ((Thl»  iRpllPi 


Sit  loriaird  Nili  1. )) 

S|  TRAVELS  THROUSH  THE  BUHPER  AT  lil  I,  <  V. 


CONDITIONS  AFTER  IMPACT 


FIGURE  3  Assumed  Penetrating  Geometry 
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As  the  impact  shocks  travel  from  their  origin  (i.e.  the  point 
of  impact),  they  will  in  general  be  diffracted.  However,  for  very  thin 
bumper  plates  we  neglect  diffraction  effects  and  ass\ame  a  one -dimensional 
wave  model  as  shown  in  Figure  3,  (b)  and  (c).  In  this  case  it  is  implied 
that  the  bumper  provides  the  necessary  retaining  walls  to  permit  the  as¬ 
sumption  of  one -dimensional  flow.  Furthermore,  shearing  stresses  are  ne¬ 
glected. 

When  the  bimiper  shock  reaches  the  far  surface  of  the  bumper 
plate,  it  is  reflected  as  a  rarefaction  ^  which  proceeds  back  through  the 
bximper-pellet  flow  at  a  velocity  (u  -  ap,  ag)  where  u  is  the  particle  ve¬ 
locity  to  the  right,  again  in  space-fixed  coordinates.  ((This  rarefac¬ 
tion  g  will  move  to  the  left  or  to  the  right  in  space-fixed  coordinates 
according  as  7g  is  greater  than  or  less  than  2--See  Section  VI.))  JR  will 
eventually  overtake  and  cancel _^p,  thus  if  the  pellet  is  long  enough  rela¬ 
tive  to  the  bumper  thickness,  a  portion  of  the  pellet  would  remain  intact 
after  penetrating  the  bxmiper.  If  ^  actually  moves  to  the  right  in  space- 
fixed  coordinates  (requiring  that  73  <  2)  then  the  compressed  regions  of 
the  pellet  and  bumper  can  protrude  through  the  bumper;  in  this  case  the 
restraint  afforded  by  the  bumper  will  cease  to  exist  at  the  rear  face  of 
the  bumper  and  a  cylindrical  rarefaction  -l-R  will  be  formed  as  this  re¬ 
straint  is  relaxed.  The  resulting  model  is  shovm  in  Figure  9{'b);  both 
the  rarefaction  R  and  4- R  accelerate  the  material  to  escape  speed.  In 
this  case  the  inward  facing  rarefaction  waves  4 R^  are  accelerating  par¬ 
ticle  flow  radially  and  attenuating  the  pellet  shock  Sp;  these  cylindri¬ 
cal  wave  fronts  proceed  radially  inward  at  the  local  velocity  of  sound, 
ap,  ag,  and  are  assumed  to  cancel  in  all  regions  where  they  intersect 
Thus  R  accelerates  the  flow  axially,  while  4  R^,  produces  radial  accelera¬ 
tion;  the  amount  and  nature  of  material  accelerated  in  both  directions  de¬ 
pends  on  the  relative  velocities  and  geometries  of  these  two  fronts. 

Thus  it  is  seen  in  this  assumed  model  that  on  impact  the  kinetic 
energy  of  the  pellet  is  transferred  partially  into  kinetic  energy  of  the 
affected  regions  of  the  bumper  plate  and  partially  into  thermodynamic  en¬ 
ergy  of  both  pellet  and  bumper  material.  The  energy  transfer  is  effected 
through  the  transmitted  and  reflected  shock  waves.  A  large  portion  of 
the  thermodynamic  energy  imparted  to  the  pellet  and  displaced  bumper  mate¬ 
rial  is  ultimately  transferred  to  kinetic  energy  through  a  series  of  ex¬ 
pansion  waves  originating  at  the  boundaries.  The  remainder  of  the  thermo¬ 
dynamic  energy  is  involved  with  the  vaporization  energies  of  the  materials. 

((The  case  of  7g  >  2  in  which  the  rarefaction  wave  accelerates 
the  material  to  less  than  escape  velocity  and  subsequent  expansion  occurs 
at  the  rear  surface  of  the  bumper  where  the  cylindrical  constraint  is  re¬ 
leased  requires  a  somewhat  more  complicated  model.  This  case  is  shown  in 
Figure  8 . ) ) 


3rthur  Sl.Hittle.Snr. 
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III.  THE  CONSERVATION  RELATIONS 


The  analysis  is  based  on  negligible  shear  forces,  i.e.  a  com¬ 
pressible  fluid  is  assumed.  The  mass  and  momentum  conservation  relations 
may  be  written  in  the  usual  form;  where  it  is  assumed  that  the  fluid  re¬ 
mains  ideal: 


Pb/b  =  ^B^^B  -  ^B^ 


^B  ^  Pb  '^B^  =  ^B  ^  Pb^^B  -  '^b)^  ') 
o  o 


across  the  bumper  shock  (5) 


and 


Pp  ’'^P  =  ^P^'^P  -  '"P^ 

o 


Pp  +  Pp  Wp  Pp  +  ■  ^p'^  ; 

o  o 


across  the  pellet  shock  (^) 


The  symbols  are  self-evident  from  Figure  . 


Figure  4 
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Equations  5  and  4  can  be  combined  and  algebraically  arranged  to 
yield  the  normal  expressions  for  shock  wave  and  particle  velocity; 


The  pressure  relations  may  be  derived  directly  also: 


^B  '^B  ■  ^B^'^B  "  ’^B^ 

o  o 


=  ^B^B^^B  -  =  Pb 


and  since  the  undisturbed  pressures  may  in  all  cases  be  regarded  as  negli¬ 
gible  in  comparison  with  the  shocked  pressures,  using  (2)  we  may  write 
simply: 


P  =  P  = 
P  B 


0-,  W  u„ 
B  B  B 


^B  Vp 

o 


(7) 


By  using  (5)  and  (ii)  we  may  write  the  appropriate  conservation  of  energy 
relations  as; 


e 


B 


1/2 


(8) 


e 


P 


1/2  \  Pp  + 


(9) 


where  the  latent  heat  of  vaporization  has  been  neglected.  ((This  neglect 
of  the  heat  of  vaporization  implies  that  we  are  considering  very  high 
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energy  processes  such  as  are  experienced  in  meteoric  impact;  the  experi¬ 
mental  data  obtained  in  hypervelocity  firings  pertain  to  lower  velocity 
impacts  in  which  this  simplification  is  inappropriate . ) )  To  sufficient 
accuracy  these  may  be  written  as: 


(10) 


(11) 


For  polytropic  gases,  the  enthalpy  may  be  written  as  ( v/y-l ) • (p/p )  and 
the  standard  Rankine-Hugoniot  pressure -density  relations  derived: 


rp  +  1 


7b  +  1 


7p  -  1 


7t,  +  1  P 
'P  o 


7p  -  1  Dp 


-  1 


’'b  -  ^ 

o 


’'b  "  ^  ^o 


(12) 


1  D. 


-  1 


which  defines  the  limiting  density  ratio  for  very  strong  shock  waves  as: 


lim 


oo 


0p\ 


7p  +  1 


lim 


QD 


"B 


’'b  -  ^ 


(15) 


Based  on  the  ass\miption  that  the  impact  shocks  are  very  strong,  (l5)  may 
be  applied  to  yield  the  following  relations: 

a.  Particle  velocity  ratio: 


/  ^P  "  " 


(14) 
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"b.  Particle  velocity  ‘behind  the  shock:  (combine  Equations  1  and  l4.) 


1  + 


^  "  ^o 
7p  +  1  Pp 


(15) 


c.  Shock  wave  velocities  (combine  Equations  5  and  15)' 


7b  +  1 


7p  +  1 
2 


(16) 


d.  Pressure  in  shocked  states: 


7b  +  1 

^P  "  ^B  "  2  ^B^  “ 


1  + 


7t3  +  1  B 
_B _ o 

7p  +  1  Pp 


(17) 


e.  Temperature  in  shocked  states: 

This  may  be  obtained  by  combining  the  strong  shock  relations 
with  (10)  and  (11),  using  (17),  and  assuming  the  ideal  gas  relation 
e  =  c  T: 

V 


T„  -  T. 


’'b  -  ^ 


B  Bo  2Rb 


1  + 


7p  +  1  Pp 


or  to  sufficient  approximation: 


’'b  -  ^ 


"B  2R, 


1  + 


^B_ii  % 

7p  +  1  Pp 


(18) 
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and  similarly: 


f .  The  velocity  of  sound  in  the  shocked  states  may  be  derived  from  the 
polytropic  relation  a^  =  7p/P»  yielding: 


Equations  li<-  to  22,  inclusive ,  completely  define  the  impact  induced  thermo¬ 
dynamic  states . 
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IV.  SOME  NUMERICAL  EXAMPLES 


In  the  subsequent  analysis  of  three-dimensional  motion,  it  is 
of  prime  importance  to  .  know  the  practical  range  of  conditions  to  be  en¬ 
countered  behind  the  impact  shocks.  In  the  condensed  states,  such  as  are 
encountered  behind  the  impact  shocks,  considerable  uncertainty  arises  as 
to  the  value  of  7  to  be  chosen.  In  the  subsequent  calculations,  quanti¬ 
ties  will  be  expressed  in  terms  of  7  under  the  assumption  that  7p~ 
so  that 


rp  +  1 


1, 


but  since,  y^,  y„  may  approach  unity,  the  inequality 
r  ii 


will  be  observed. 

Various  quantities  estimated  in  Equations  li^  to  22  have  been 
evaluated;  these  are  plotted  in  Figures  6,  and  7* 


3ttluir  Sl.'Uittlc.Urtr. 
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IMPACT  VELOCITY,  km/sec 

FIGURE  5  Impact  Pressures  from  Equation 
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V.  DISCUSSION  OF  RESULTS 


One  of  the  most  important  assumptions  (aside  from  that  of  ideal 
fluid  behavior)  is  the  use  of  the  strong  shock  approximation,  i.e.  the 
use  of  the  limiting  density  ratio.  The  rigorous  justification  of  such  an 
approximation  is  obviously  not  possible .  Its  effect  must  be  estimated  in 
relation  to  each  parameter. 


A.  PRESSURE  RESULTS 


The  pressure  is  plotted  in  Figure  5  for  a  range  of  cases;  the 
effect  of  initial  density  ratio  may  be  seen  by  comparison  of  the  magne¬ 
sium  into  aluminum  case  (Mgp/Alg)  with  the  copper  into  aluminum  case 
(Cup/Alg);  one  finds  a  doubling  of  pressure  ratio  in  the  case  of  the  heav¬ 
ier  metal  impacting  pellet. 

In  these  plots  y  has  been  used  as  a  parameter  multiplying  pres¬ 
sure  .  It  may  be  observed  that  for  condensed  media,  the  appropriate  value 
of  y  is  generally  accepted  as  3  (see  references  U  and  8).  An  interesting 
comment  may  be  made  on  this  value  by  examining  Equation  17  for  the  impact 
of  similar  materials,  in  which  case 


P 


7  +  1 
2 


2 


which  with  7=5  reduces  to  the  well-known  expression  used  by  many  authors 
for  Impact  pressure: 


P 


It  may  be  seen  by  consideration  of  the  exact  Rankine -Hugoniot  equations 
relating  pressure  ratio  to  density  ratio  across  a  shock,  that  for  large 
values  of  y  the  limiting  density  ratio  is  rapidly  attained.  Thus  for 
7=3  one  would  expect  reasonable  pressure  predictions  using  the  strong 
shock  assumption.  Bjork  (Reference  3,  Figure  4)  quotes  a  computed  pres¬ 
sure  variation  for  iron  impacting  onto  iron  as  p  =  0.120  where  p  is 

in  megabars  and  V  in  km/sec .  This  value  has  been  included  in  Figure  5 
for  comparison  with  the  Fep/Fe^  pressures  predicted  by  Equation  17 .  In 
plotting  Bjork' s  formula,  the  pressures  have  been  reduced  to  correspond 
to  the  case  7=3.  The  agreement  with  the  predictions  of  Equation  I7  may 
be  observed  to  be  good. 

Experimental  measurements  of  shocked  pressures  in  metals  are  few.  Among 
the  most  extensive  studies  were  those  of  the  Los  Alamos  Laboratory 
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(References  5>6,7)"  Unfortunately,  driver  free -surface  velocity  is  not 
given  in  their  reports,  so  that  only  qualitative  comparison  can  he  made. 
In  this  case  it  can  be  seen  that  based  on  a  7  =  3,  and  assuming  7p  =  /g, 
an  impacting  steel  plate  will  generate  1  megabar  shocks  in  a  brass  speci¬ 
men  at  an  impact  velocity  of  approximately  1^,000  fps  and  I.5  megabar 
shocks  at  an  impact  velocity  of  approximately  17,000  fps. 


B.  TEMPERATURE  IN  SHOCKED  STATES 

Temperatures  in  the  shocked  state  as  given  by  Equations  18  and 
19  have  been  plotted  in  Figure  6.  It  is  interesting  to  note  the  simple 
form  of  the  temperature  relation  when  similar  materials  collide; 


’’b  =’'p 


The  temperatures  calculated  directly  from  the  analytical  relations  do  not 
take  into  account  the  energy  absorbed  in  breaking  the  crystal  lattice  and 
vaporization.  Using  the  values  quoted  by  the  Russian  workers  (Reference  9^ 
page  307 ) : 


Medium 

Energy  Density  of  Crystal  Lattice 
Including  Latent  Heat  of  Fusion 

Vaporization 
Ene rgy 

ergs 'gm 

ergs/gm 

Sand 

5  X  10^ 

Granite 

7  X  10^ 

2  X  10^° 

Aluminum 

X  10^ 

lO^l 

Iron 

3  X  10^  . 

7  X  10^° 

The  predicted  temperatures  may  be  reduced  to  allow  for  the  vaporization 
energy,  where  the  decrease  in  the  latent  energies  with  temperature  must 
be  properly  taken  into  account.  When  these  calculations  are  performed, 
it  is  found  that  negligible  decreases  in  temperature  occur  at  the  higher 
temperatures  where  this  theory  is  essentially  applicable;  the  main  modi¬ 
fications  occur  at  the  lower  end  of  the  velocity  range,  and  here  the 
Deybe  model  along  with  the  Mie -Gruhlesen  equation  of  state  as  applied  in 
References  5^6,7,  yields  temperatures  approximately  one  order  of  magni¬ 
tude  lower  than  predicted  by  the  theory  presented  here.  Generally  speak¬ 
ing,  it  would  be  expected  that  above  10  km/sec  the  vaporization  energy 
may  be  neglected  and  the  predictions  of  the  simple  theory  considered 
valid  other  than  as  limited  by  the  ideal  type  behavior  assumed. 

The  simple  expression  quoted  above  for  impact  temperatures  when 
similar  materials  collide,  shows  that  the  temperature  rise  (7  assimied  the 
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same  in  all  cases)  varies  between  materials  directly  as  the  molecular 
weight.  Thus  the  three  curares  A1/A1,  Fe/Fe,  and  W/W  of  Figure  6  reflect 
simply  the  molecular  weight  changes.  The  more  complex  expressions  in¬ 
volved  when  different  metals  collide  yields  the  simple  relation  between 
pellet  and  bumper  temperatures: 


(Mol.Wt.)j, 

(Mol.Wt.)_  '^P 
o'  o 


The  extremely  high  temperatures  predicted  naturally  indicate 
the  break-down  of  the  simple  ideal  theory  used  here .  First  as  the  high 
temperatures  are  approached,  electrons  will  be  stripped  off  the  atoms, 
yielding,  of  course,  a  lower  gas-kinetic  temperature,  and  requiring  the 
application  of  a  Thomas-Fermi  equation  for  the  condensed  state.  In  the 
second  Instant,  with  the  high  temperatures  achieved,  it  is  apparent  that 
radiative  processes  will  tend  to  dominate  the  early  phases.  ( (Bremsstrah- 
lung  and  absorption-re -emission  may  well  prove  to  be  important  mechanisms 
for  transporting  energy  from  the  shocked  interior  to  the  radiating  sur¬ 
face.))  The  gas  flow  processes  will  occur  relatively  at  a  much  slower  rate 
until  the  extreme  temperatures  have  been  dissipated.  Both  the  above  com¬ 
ments  lead  to  the  added  complexity  of  non-negligible  departures  from 
local  thermodynamic  equilibrium,  while  inherent  in  the  simple  model  here 
is  the  assumption  of  a  Maxwellian  distribution. 

Experimental  evidence  to  date  on  impact  flash  temperature  is 
practically  non-exlstant;  a  few  spectra  have  been  taken,  as  well  as  some 
long  wave  length  radiometric  observations,  but  only  general  qualitative 
remarks  can  be  made  on  these  observations.  The  intensity  of  the  radia¬ 
tion  throughout  the  visible  and  near  visible  regions,  even  at  low  impact 
velocities  (6-7  km/sec),  is  attested  to  by  the  clarity  of  the  image  left 
on  an  exposed  photographic  plate  (shutter  left  open),  the  motion  of  the 
pellet  being  frozen  at  the  impact  condition  by  its  own  impact  flash. 

Rather  intense  X-band  pulses  have  been  recorded  on  impact  in  this  velocity 
range;  a  black-body  calculation  would  indicate  a  radiant  intensity  (based 
on  a  pellet  length  of  shocked  gas  ^  5)  of  the  order  of  a  few  microwatts 
in  this  frequency  range  for  temperatures  of  the  order  of  lO^^K.  Observa¬ 
tions  have  shown  considerably  larger  emissions  az  these  long  W8.ve  lengths 
(3  cm),  and  these  have  been  attributed  tentatively  to  Bremsstrahlung.  To¬ 
tal  thermal  intensities  may  be  shown  to  be  quite  large  based  on  the  ideal 
model  (of  the  order  of  10®  to  10^2  watts).  Thus  the  tentative  interest¬ 
ing  observation  may  be  made  that  as  real  gas  effects  are  introduced  and 
the  gas-kinetic  temperature  drops,  the  power  emitted  tends  to  Increase  in 
the  long  wave  length  regions.  At  the  same  time,  of  course,  departures 
from  local  thermodynamic  equilibrium  become  important. 
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VI.  PELLET  DYNAMICS  AND  WAVE  GEOMETRY 


The  one -dimensional  status  may  be  assumed  to  exist  •until  the 
bumper  shock  has  reached  the  far  side  of  the  bumper  plate,  at  which  time 
it  will  be  reflected  as  a  rarefaction  centered  about  this  rear  face.  Un¬ 
der  the  assijmptlons  of  the  theory,  the  rear  face  would  then  be  accelerated 
instantaneously  to  escape  velocity,  and  the  rarefaction  spreading  out  as 
it  progresses  would  travel  back  through  the  bumper-pellet  at  a  velocity 
'^B,P  “  Two  conditions  are  of  interest  here.  First,  if  the  front  of 
the  rarefaction  travels  upstream  through  the  condensed  media  at  a  velocity 
greater  than  the  local  particle  velocity,  then  the  rarefaction  will  never 
pass  downstream  of  the  rear  bumper  surface,  and  in  effect  we  would  have 
the  media  issuing  from  a  cylindrical  orifice  with  the  initial  front  ve¬ 
locity  approaching  escap)e  speed  and  gradients  appropriate  to  the  centered 
rarefaction,  with  the  gas  column  expanding  laterally  to  adjust  to  the 
boundary  conditions  of  free  space.  This  condition  is  illustrated  crudely 
in  Figure  8. 


FIGURE  8 


The  analytical  condition  for  which  this  state  of  affairs  may  be  expected 
can  be  expressed  by  comparing  Equations  15  and  20,  and  it  may  be  seen 
that  the  requirement  is: 
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1) 

—  >  1 


or: 


7b  >2 


If  7b  >  2,  then  it  is  apparent  that  the  front  of  the  rarefac¬ 
tion  wave  Will  be  swept  past  the  rear  face  of  the  bumper^nd  th?  restraint 
on  the  boundary  of  the  shocked  volume  contained  between  Sg  and  suddenly 
released.  Under  these  conditions,  it  would  appear  a  reasonable  first -try 
approximation  to  assume  that  the  bianper  plate  is  removed  instantaneously 
from  the  picture,  having  established  in  the  pellet  a  pulse  of  one -dimensional, 
shocked,  condensed  media,  of  width  somewhat  less  than  twice  the  bumper 
plate  thickness.  Thus  we  imagine  that  at  the  instant  the  bumper  impact 
shock  reaches  the  far  surface  of  the  bumper  plate,  all  restraint  is  removed 
on  the  pellet  boundaries,  and  the  wave  geometry  appropriate  to  this  assump¬ 
tion  is  illus^ated  in  Figure  9(a).  At  this  instant  a  plane  centered  rare¬ 
faction  wave  starts  upstream  through  the  pellet  and  must  overtake  the 
pellet  shock  §b  eventually  (if  the  pellet  is  sufficiently  long). 

In  addition,  a  centered  cylindrical  rarefaction  *  R^  must  start 
at  the  free  boundary  and  progress  inwards  tovrards  the  center  of  the  pellet, 
imparting  radial  acceleration  to  the  condensed  media.  After  some  time,  the 
wave  system  in  the  pellet  will  have  assumed  the  form  shown  in  Figure  9Cb) • 

Radial  outflow  will  be  initiated  by  the  passage  of  the  impact 
shock  and  stopped  by  the  arrival  of  the  end  rarefaction  wave.  Expansion 
waves  from  the  cylindrical  boundaries  will  be  propagated  throughout  the 
field,  overtaking  and  decaying  at  least  part  and  possibly  all  of  the  im¬ 
pact  shock  before  the  arrival  of  the  end  rarefaction.  This  diffraction  at 
the  edge  of  the  impact  shock  is  illustrated  in  Figure  9(b)*  The  pellet 
flow  regions  may  be  divided  then  as  follows: 

a.  A  one -dimensional  field  of  axial  high  velocity  gas  flow 
behind  the  plane  portion  of  the  end  rarefaction  (Region  1 
in  Figure  9(b)). 

b.  A  one -dimensional  field  of  axial  lower  velocity  flow  behind 
the  planar  portion  of  the  impact  shock  and  the  boundaries 
of  •IRc,  (Region  2  of  Figure  9(b)). 

c.  A  more  or  less  radial  outflow  field  ahead  of  the  front  R-^^, 
established  by  J*Rc  (Region  5  of  Figure  9(b)). 
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d.  An  interaction  field  Joining  the  sixial  flow  region  (l)  to 

the  radial  flow  r^ion  (5),  established  by  the  interaction 
between  and  R^. 

e.  An  interaction  field  Joining  the  radial  field  (5)  with  the 
axial  field  (2)  and  established  by  the  interaction  of  edge- 
rarefactions  with  the  impact  shock. 


The  detailed  solution  of  each  of  these  fields  is  beyond  the 
scope  of  this  present  note  and  will  be  treated  in  subsequent  reports. 
There  are,  however,  some  interesting  simple  points  that  might  be  made 
without  too  much  complexity  of  analysis.  In  the  first  instance,  one  is 
led  to  wonder  how  soon  the  impact  shock  will  be  decayed,  either  by  cylin¬ 
drical  edge  disturbances  or  by  the  overtaking  of 

The  edge  disturbance  may  be  considered  in  a  simple  geometrical 
manner  with  reference  to  Figure  10,  where  the  conditions  are  shown  an  in¬ 
finitesimal  time  after  the  release  of  side-wall  restraint.  An  axis  sys¬ 
tem  is  shown  fixed  in  the  piellet  at  the  time  of  release. 
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The  expansion  disturbance  originating  at  (O,  rg)  will  be  translated  behind 
the  shock  at  velocity  Up,  and  will  penetrate  in  and  Interact  with  the 
shock  to  a  depth  dy.  The  differential  equation  defining  this  inner  point 
of  interaction  is: 


dx  Wp 


which  can  be  integrated  (assxaning  ap,  Wp,  Up  constant)  to  yield: 


which  on  substitution  of  Equations  l4,  15,  I6,  and  21  yields: 


y 

X 


These  sidewaves  wilj^have  penetrated  to  the  pellet  axis  (and  hence  attenua¬ 
ted  the  whole  wave  Sp)  at: 


which  for  7p  =  3,  is  approximately  1.5  radii.  Thus  we  are  in  effect  saying 
that  the  shock  diffraction  effects  are  confined  to  a  small  portion  of  the 
planar  front  for  bodies  whose  lengths  are  less  than  their  diameters. 

The  end  originating  rarefaction  R  will  overtake  and  interact 
(tending  to  cancel)  the  pellet  impact  shock  wave.  This  may  be  expected  to 
be  the  most  important  cancellation  process,  where  lateral  dimensions  are 
of  at  least  the  same  order  of  magnitude  as  longitudinal  pellet  dimensions. 
The  position  in  the  pellet  for  which  overtaking  will  occur  may  be  computed 
from  the  simple  geometry  of  Figure  11 . 


J 

1 

I 

I 

I 

I 

I 

I 
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R 


Figure  11 


In  Figure  11,  denotes  the  length  of  pellet  which  must  enter 
the  shock  wave  Sp  from  the  time  the  rarefaction  wave  is  first  formed  until 
it  overtakes  the  shock.  The  pellet  of  length  will  he  compressed  to 
Xi(l  -  Up/wp)  on  passing  through  the  shock  wave.  Equating  the  times  for 
the  shock  and  the  rarefaction  to  merge,  one  has 


-1  =  2_  (1  _  P  (1  .  ^)  +  Jl.  (1  _ 

^  ^P  ^P  '^P 


Upon  substitution  of  the  appropriate  relations  from  Equations  l4  to  22, 
one  obtains 
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/Pg  (rp  +  1) 

o _ 

o 


-  1) 

V^'P  ■ 


+  1 


2r. 


rp  -  1 


-  1 


wherein  substitution  of  7  =  7ri  =  7  yields 

r  B 


(23) 


This  function  has  been  plotted  in  Figure  12  as  a  function  of  7 
and  initial  density  ratio,  and  it  may  be  observed  that  overtaking  occurs 
in  quite  a  short  distance ,  probably  within  ten  times  the  bumper  plate 
thickness,  in  the  most  probable  region  of  interest.  The  interaction  of 
an  overtaking  rarefaction  with  a  shock  wave  has  been  discussed  by  several 
authors  (Reference  11),  but  is  complicated  in  this  case  by  degeneracy 
into  elastic  shocks. 


FIGURE  12  Overtaking  Point  Measured  from  Initial  Position 
of  Sp  at  Time  Sg  Reaches  far  Side  of  Bumper 
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VII .  RADIAL  FLOW 


As  previously  stated,  the  band  of  highly  shocked  gas  between 
the  wave  fronts  and  R  will  start  to  expand  radially  as  soon  as  the 
sidewall  pressure  is  released,  i.e.  at  the  rate  the  boundary  is  swept  by 
the  rear  bumper  surface.  For  the  first  approximation  we  have  assumed  in¬ 
stantaneous  release  of  the  restraint  on  the  sidewall  so  that  radial  out¬ 
flow  starts  at  the  shock  wave  «Sp.  The  detailed  geometry  was  shown  in 
Figure  9(13).  The  reference  axis  chosen  is  fixed  in  the  pellet  at  the  un¬ 
disturbed  end.  With  respect  to  this  axis  system,  ^Sp  induces  a  particle 
velocity  up  in  the  negative  direction.  The  inward  traveling  rarefaction 
nI  Rj,  accelerates  the  flow  approximately  radially,  and  since  it  is  cen¬ 
tered  with  its  origin  at  r  =  r^,  then: 
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where  it  is  assumed  that  the  polytropic  exponent  y  may  be  approximated  by 
a  constant  average  value.  In  practice  it  is  known  that  condensed -media 
on  expansion  must  follow  a  more  complicated  form  where  y  =  7(p,T),  y  as¬ 
suming  its  Isentroplc  value  only  after  considerable  expansion  (Stanyuko- 
vich  suggests  for  p  <  10^  atmospheres) .  This  problem  will  be  considered 
in  detail  in  a  subsequent  paper. 

Since  the  flow  in  this  region  is  translating  with  a  velocity  u  , 
the  envelope  of  the  expansion  products  will  be  a  cone  with  semi-vertex 
angle : 
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(the  subscripts  defining  the  envelopes  for  the  two  different  materials). 
See  Note  2  of  foreword. 

Using  the  previously  derived  relations,  this  may  be  written  as: 
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for  the  bumper  material,  and; 


29 


tan 


-1 


ilp 

7_  -  1  +  1 

'P  'P 

o 


(25) 


for  the  pellet . 

Assuming  that  for  the  very  high  pressures  and  temperatures  in¬ 
volved  behind  the  impact  shocks,  metal  exponents  are  approximately  identi¬ 
cal,  i.e.  7p  =  7-q)  these  functions  have  been  plotted  in  Figure  15 .  It  is 
interesting  to  note  that  Equations  24  and  25  permit  determination  of  mean 
7's  from  spray  angle  envelope  measurements  (assuming  the  expansion  law  is 
known).  It  is  to  be  rr :ed  that  the  spray  angle  is  quite  sensitive  to 
both  polj-tropic  gas  exx>onent  and  the  combination  of  impacting -bumper  mate¬ 
rials;  a  light  pellet  impacting  on  a  dense  bumper  (magnesium  onto  stain¬ 
less  steel)  doubles  the  spray  angle  for  the  pellet  from  the  reverse  condi¬ 
tion  (copper  pellet  onto  an  aluminum  bumper). 

It  may  be  observed  that  the  values  yielded  by  Equations  23  and 
24  appear  to  be  high  when  compared  to  existing  experimental  data  from  hy¬ 
personic  range  firings,  e^.'en  assuming  the  large  \'alues  of  7.  Two  comments 
may  be  made  here,  first  the  simple  theory  advanced  here  does  not  account 
for  vaporization  energy  and  condensation  effects  on  undergoing  the  sudden 
type  expansion  associated  with  the  centered  wave.  In  the  second  instance, 
experimental  results  to  date  have  been  at  too  low  velocity  to  provide  any 
cross  checking;  impact  energies  were  not  sufficiently  high  to  meet  the  as¬ 
sumptions  inherent  in  this  simple  model,  i.e.  vaporization  energy/shock 
enthalpy  much  smaller  than  unity. 

Assigning  cylindrical  flow,  (such  as  would  predominate  in  Region 
3  of  Figure  9(4i))  the  eqiiations  to  be  solved  are: 
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7  =  7(P>T) 


with  boundary  conditions  those  of  a  centered  wave  at  r  =  r  ,  radial  flow 
durations  being  bounded  by  the  interaction  of  the  rarefact?ons .  These 
equations  along  with  the  appropriate  energy  relation  can  be  programmed 
for  numerical  solution,  including  with  somewhat  considerable  increase  in 
complexity,  detailed  analysis  of  the  interaction  regimes.  If  a  constant 
mean  value  of  7  can  be  assumed  over  some  expansion  region,  then  the 
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classical  self-similar  solutions  (Sedov  11,  Stanyukovich  9)  can  be  ap¬ 
plied  directly  to  give  radial  distribution  of  velocity,  density,  and 
energy. 
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VIII.  Discussion 


The  present  note  Is  intended  to  suaeserite  aooe  Inltisl  thlR?slnif 
on  an  analj-tical  approach  to  the  break-up  of  a  meteor  on  Impact  with  s 
thin  flat  plate.  Periodic  extensions  a.nd  oodiflcations  of  th*  approach 
will  he  made  in  light  of  more  extensive  work;  Indeed,  •.  •.♦•n  at  the  tla>-  of 
issuance  some  of  the  work  has  been  advanced  well  beyond  the  sts^s  de¬ 
scribed  herein.  However,  it  is  felt  worthwhile  to  circulate  the  present 
write-up,  since  it  does  contain  consldertible  detail  that  will  not  he  re¬ 
peated  in  future  publications. 

In  this  analysis  three  main  phases  of  the  protlcs  haw  been  dis¬ 
cussed,  some  in  considerably  more  depth  than  oth<*ra.  One  1*  naturally 
led  to  wonder  in  summary,  the  states  of  e&^h  phase. 

In  relation  to  the  first  phase,  'h  predictions  of  pressure  and 
temperature  in  the  Impact  shock  states  would  nrpeur  to  te  satisfactory  as 
a  first  approximation,  barring,  of  course,  significant  electron  stripping 
effects.  The  comparison  with  BJork's  computed  curw  indicates  *0!te  quite 
reasonable  agreement  with  more  exact  computations.  These  ofc»erv?»tlons 
are  all  based  on  the  use  of  the  condensed  a'^dia  •>  ■  5.  ’to  co«f*rSsos  can 
be  made  with  any  experimental  or  predicted  temper'tturc  effects;  ‘he  los 
Alamos  work  considered  only  the  solid  state  propagation  of  shocks  and 
falls  well  below  the  applicability  of  this  theory.  Electron  stripping  will 
lower  the  temperature  and  raise  the  pressure,  and  this  is  presently  being 
included  in  estimation  of  the  impact  states. 

In  relation  to  the  second  phase,  considerably  s»re  effort  «ust 
be  expended  on  determining  the  radiant  energ;.’  fl'ux  and  its  effect  on 
shock  induced  and  expansion  states.  Cr-de  estimates  wopld  indicate  ra¬ 
diative  losses  throughout  the  spectr  ur  in  t?rc  10'  to  ID"*  watts  per  square 
centimeter  range.  These  may  be  expected  to  lower  appreciably  the  therso- 
dynamic  state  at  which  expansion  takes  place  as  well  as  introducing  dejufcr- 
tures  from  local  thermodynamic  equilibrium.  This  chase  is  perhaps  in  the 
least  satisfactory  condition  of  all. 

Given  the  thermodynamic  states,  the  expansion  flo'w  characteris¬ 
tics  can  be  computed  readily.  Simple  boundary  controlled  conditions  have 
been  presented  in  this  note,  ur.i  the  more  detailed  solutions  to  tix*  expan¬ 
sion  flows  will  be  presented  in  another  note  soon  to  be  published.  The 
more  complex  problem  of  expansion  from  the  condensed  state  (with  a  >'  ■  y(p»T)) 
remains . 


Experimentally,  shocks  have  been  propagated  through  columns  of 
gas,  and  the  diffraction  and  outflow  characteristics  obserwd.  These  stud¬ 
ies  tend  to  establish  the  validity  of  the  basic  model  and  provide  some 
checks,  crude  though  they  may  be,  on  the  boundary  flows.  Future  experi¬ 
mental  work  will  include  the  measurement  of  impact  radiation  and  spray 
characteristics  in  the  McGill  hypersonic  firing  range. 
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